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Many transition metal compounds show a saturation of the electrical resistivity at high temper- 
atures, T, while the alkali-doped fuUerenes and the high-Tc cuprates are usually considered to show 
no saturation. We present a model of transition metal compounds, which shows saturation, and a 
model of alkali-doped fuUerenes, which shows no saturation. The electron scattering is assumed to 
be due to interaction with phonons. The properties of these models are determined by performing 
quantum Monte-Carlo calculations. To analyze the results, as well as earlier results for the high-Tc 
cuprates, we use the f-sum rule. We demonstrate that the f-sum rule leads to a natural upper limit 
for the resistivity at large T. For some systems and at low T, the resistivity increases so rapidly 
that this upper limit is approached for experimentally accessible temperatures. The resistivity then 
saturates. For a model of transition metal compounds with weakly interacting electrons, the upper 
limit corresponds to an apparent mean free path consistent with the loffe-Regel condition. For a 
model of the high Tc cuprates with strongly interacting electrons, however, the upper limit is much 
larger than the loffe-Regel condition suggests. This upper limit is not exceeded by experimental 
resistivities. The experimental data for the cuprates are therefore consistent with saturation. After 
saturation the resistivity normally grows slowly. The alkali-doped fuUerenes can be considered as 
systems where saturation has happened already for T — 0, due to orientational disorder. We show, 
however, that for these systems the resistivity grows so rapidly after "saturation" that this concept 
is meaningless. This is due both to the small band width and to the coupling to the level energies 
of the important (intramolecular) phonons in the fuUerenes. 



I. INTRODUCTION 

The electrical resistivity of metals is often described 
in a semiclassical picture, where an electron on the av- 
erage travels a mean free path I before it is scattered by 
a phonon, an impurity or another electron. Assuming a 
spherical Fermi surface, the resistivity p can be expressed 
in terms of / as 



(1) 



where kp is the Fermi wave vector. Alternatively, if we 
know the resistivity experimentally, we can deduce an ap- 
parent mean free path from Eq. (1). For a good metal, 
I is typically several hundred A or more. As the temper- 
ature T is increased, p increases. Normally, it is found 
that p{T) ^ T for T larger than some fraction of a typical 
phonon energy. This is due to the increased scattering by 
phonons, and it corresponds to a reduction of I. Never- 
theless, at the melting point, I is still typically very much 
larger than the separation d of two neighboring atoms. 
An example of this behavior is given by Cu in Fig. 1. 

In the 1970's a number of exceptions to this behavior 
were found. ^ In particular for several A15 compounds, 
such as NbaSb and NbaSn, it was found that p increases 
very rapidly with T for small T, leading to very large 
values already for temperatures of the order of a few hun- 
dred K. At these values of T, the slope of p{T) is strongly 
reduced. This is shown in Fig. 1, where the resistivi- 
ties of NbaSb and Cu are compared. This was described 
as "resistivity saturation".^ Interestingly, it was found 
that saturation happened when I ~ d, the loffe-Regel 
condition.^ The corresponding resistivity is also shown in 



Fig. 1. During the 1970's and early 1980's many exam- 
ples of this were studied, and saturation of the resistivity 
when / ~ d was considered a universal behavior. 
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FIG. 1. Resistivities of Cu and NbaSb 
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^ The figure also 
saturation resistivity for NbsSb, ob- 



shows the loffe-Regel 
tained by assuming that the mean free path I in Eq. (1) is 
equal to the the distance between the scattering centers. The 
figure illustrates that for NbsSb the resistivity saturates at 
roughly the value expected from the loffe-Regel criterion.^ 



In a semiclassical picture, this behavior may be ex- 
pected. It may seem that the worst that could happen 
is that an electron is scattered at every atom. We would 
then expect / d to be fulfilled. This argument is, how- 
ever, not convincing. In the semiclassical theory, it is 
assumed that an electron travels through the solid with 



1 



a well-defined k-vector between the scattering events. If, 
however, I ^ d, it is not possible to define k, and the 
theory breaks down.^ A proper theory of saturation is 
therefore needed. A number of theories have been put 
forward, but no theory has been generally accepted. 
Due to the break-down of the semiclassical theory when 
I ~ d, the concept of a mean free path itself becomes 
questionable for such small values of Z. In this case we 
use Eq. (1) as a definition of the (apparent) mean free 
path. 
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FIG. 2. Resistivities of Lai.93Sr().07CuO4* and RbaCeo^ 
and the corresponding loffe-Regel^ saturation resistivities.^^ 
The figure illustrates that the resistivity of these systems be- 
comes much larger than predicted by the loffe-Regel condi- 
tion. 

More recently, several apparent exceptions to resistiv- 
ity saturation have been found. In particular, this is 
the case for some strongly correlated systems, for in- 
stance the high- Tc cuprates,* and for the alkali-doped 
fullerenes.^'^° This is illustrated in Fig 2, where we show 
the resistivities of Lai.93Sro.o7Cu04 and RbaCeo together 
with the loffe-Regel resistivities.^^ Different experiments 
for alkali-doped Ceo compounds show substantial differ- 
ences, but this is not essential for the present discussion. 
The lofFe-Regcl resistivities of these two systems arc very 
large, due to the low carrier density. The figure illustrates 
that the experimental resistivities,nevertheless, greatly 
exceed the loffe-Regel resistivities. It also illustrates that 
the resistivities of these two compounds are very much 
larger than for NbsSb and other systems, which shows 
saturation according to the loffe-Regel condition. 

This shows that the semiclassical argument behind the 
loffe-Regel condition is not only questionable, but that 
it leads to wrong conclusions for the high- Tc cuprates 
and the Ceo compounds. This emphasizes the need for 
a proper theory of why saturation happens for some sys- 
tems but not for others. We also need to understand why 
saturation happens for most transition metal compounds 
when I d, although I is not a well-defined concept any 



more. 

We have earlier presented such a theory for transition 
metal compounds in a short publication,^^ and we here 
expand the arguments. We have also analyzed the rea- 
sons for the lack of saturation in the alkali-doped Ceo 
compounds, and we provide additional results here. Fi- 
nally, we have also presented results for a model of the 
high- Tc cuprates. We have therefore considered models 
of three classes of systems: i) a model of weakly corre- 
lated transition metal compounds, which shows satura- 
tion in agreement with the loffe-Regel condition, ii) a 
model of strongly correlated high- Tc compounds, which 
shows saturation but at much larger values than pre- 
dicted by the loffc-Regel conditions, and iii) a model of 
alkali-doped fuUerenes, which shows no saturation. 

We assume that in case i) and iii) the important scat- 
tering is due to the electron-phonon interaction. In a 
model Hamiltonian approach, there are two natural types 
of coupling to the phonons, either via the level energies 
(LE coupling) or via the hopping matrix integrals (HI 
coupling). In most nonionic compounds the latter effect 
should be the dominating one. As the distance between 
two neighboring atoms is changed due to the excitation 
of a phonon, the main effect should be a change of the 
hopping integrals. We study this for a model of transition 
metal compounds, referred to as the TM model. 

In molecular solids, such as the alkali-doped fullerenes, 
the situation is different. Due to the weak coupling be- 
tween the molecules, it is sensible to first calculate the 
levels of a free molecule, and then to study the weak hop- 
ping between these levels. In the alkali-doped fullerenes 
the main coupling is to intramolecular phonons. These 
phonons couple primarily to the level energies and only 
weakly to the hopping integrals between the molecules. 
We therefore study the LE coupling for a model of alkali- 
doped Cfio systems, in the following referred to as the Ceo 
model. The LE coupling may also become important for 
strongly ionic systems. 

We use a quantum Monte-Carlo (QMC) method^^ for 
calculating the current- current correlation function for 
imaginary times. A maximum entropy method^^ is then 
applied to analytically continue the response function to 
the real frequency axis. This gives the frequency de- 
pendent optical conductivity <t{u)), and the resistivity 
p = 1/(7 (ti) = 0). Since the QMC method has no sign 
problem for the models studied here, we are able to ob- 
tain rather accurate results for the resistivity. In par- 
ticular, we can establish whether or not the models we 
consider show resistivity saturation. 

To interpret the results we use a simplified approx- 
imate approach, treating the phonons (semi) classically. 
By comparing with In this method we assume that the 
phonons can be described by random static displace- 
ments of the atoms with an average amplitude that in- 
creases with T. The remaining electronic problem can 
then easily be solved quantum mechanically. This ap- 
proach is in contrast to the Boltzmann equation, where 
the electrons are treated semiclassically. The main ad- 
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vantage of this method, compared with the QMC calcula- 
tion, is that it is simple enough to allow an interpretation 
of the results. By comparing with the QMC results we 
establish the range of applicability of the semiclassical 
method for the models of interest here. 

In our semiclassical treatment, the excitation of 
phonons leads to a static variation of the level energies 
in the Ceo model and of the hopping integrals in the 
TM model. In the context of disordered systems, this is 
referred to as diagonal and off-diagonal disorder, respec- 
tively. Past work has primarily studied diagonal disorder, 
which in some respects is technically simpler. 

Traditionally, transport is described within the Boltz- 
mann theory. The Bloch-Boltzmann^^ theory starts from 
the perfectly periodic system, and treats the scattering 
mechanisms as small perturbations. This can be consid- 
ered as a theory which is valid to lowest order in 1/ {kpl)-'^ 
This further empahsizes that the Boltzmann equation be- 
comes questionable when I ~ d. Furthermore, the Zi- 
man approximation"'^* to the Boltzmann equation leads 
to p{T) ~ T for large T, i.e., there is no saturation in 
contrast to what is found experimentally for many sys- 
tems. It is then natural to look for extensions of the 
Boltzmann equations, which would extend the range of 
perturbation strengths that can be treated."'''^ We find, 
however, that in, e.g., the A15 compounds the thermally 
excited phonons even at relatively small T tend to largely 
remove the effects of periodicity. In the semiclassical 
treatment of the phonons, the momentum conservation in 
the electronic system is lost already for temperatures of 
the order of a few hundred K. We therefore consider the 
opposite limit to the Boltzmann equation, where we as- 
sume that thermal excitations have completely destroyed 
periodicity. At low T there is a Drude peak in the op- 
tical conductivity a{Lo) due to intraband transitions be- 
tween states with similar k-vectors. As T is increased, k- 
conservation is lost, and the meaning of intraband transi- 
tions is blurred, the Drude peak disappears. We therefore 
focus on the limit where there is no pronounced structure 
in cr{io) at small lj. 

We have earlier used current and charge conservation 
to obtain simple upper estimates for the resistivity of a 
metal. Here we show how the same result can be de- 
rived by using the (related) f-sum rule. This approach 
has the advantage that it can also be used to discuss 
the high- Tc cuprates,^^ and that it is convenient for dis- 
cussing the fullerenes. The approach based on the f-sum 
rule therefore provides the most convenient framework 
for analyzing the different classes of materials. 

We combine the f-sum rule with the assumption that 
the Drude peak is lost. This naturally leads to an upper 
limit for the resistivity at small or intermediate T's. If 
the initial slope of p{T) is very large, p{T) reaches this 
limiting value already for experimentally available values 
of T. At this point saturation normally happens, as is 
illustrated in our TM model. The removal of the Drude 
peak could be due to any scattering mechanism, e.g., 
electron-phonon (HI or LE coupling), electron-electron 



or disorder scattering. For the TM model considered 
here, we show in a quantum mechanical treatment that 
saturation should happen roughly when the loffe-Regel 
criterion is satisfied. This is somewhat accidental and it 
is not true for a model of the high- Tc cuprates, where 
strong correlation effects leads to a larger saturation re- 
sistivity. 

While a pronounced saturation is observed for the A15 
compounds NbsSb or NbsSn, other systems, such as Nb, 
show a weaker saturation or no saturation at all. Here we 
study a simple model of A15 compounds, referred to as 
the Nbg model, -'^ where wc include the rf-orbitals of the 
Nb atoms, put on the appropriate A15 lattice, but where 
the remaining atom (e.g., Sn in NbsSn) is neglected. This 
is compared with Nb, where the atoms are put on a bcc 
lattice. These two models then only differ with respect 
to the lattice structure. This difference leads to a smaller 
plasma frequency for NbsSb and a steeper slope of p{T). 
This leads to a much more pronounced saturation for 
NbaSb. 

Even after "saturation" has happened, p{T) tends to 
continue to grow, but at a slower rate. In this respect 
there is sometimes an essential distinction between LE 
and HI coupling. This can be best discussed using the 
f-sum rule. We show that the change of the resistiv- 
ity can be viewed as resulting from a change of the ki- 
netic energy and of the band width. These changes keep 
growing without limit with T, due to the Bose nature 
of the phonons and the lack of limitation on the number 
of phonons. The two changes work together for the LE 
coupling, but tend to compensate each other for the HI 
coupling. As a result the resistivity grows more slowly 
after "saturation" for the HI coupling and the saturation 
is more pronounced. This distinction is fairly clear-cut 
for the Ceo model. For this model, disorder leads to such 
a strong scattering, that "saturation" can be considered 
to have happened already at T = 0. Due to the LE 
coupling and the small band width, however, p{T) grows 
so rapidly after "saturation" that the concept of satura- 
tion becomes meaningless. For HI coupling, on the other 
hand, the resistivity shows a clear change in slope, even 
for the Ceo model. 

In Sec. II we present the TM and Ceo models and 
in Sec. Ill the QMC and semiclassical methods are de- 
scribed. The results are presented in Sec. IV and dis- 
cussed in Sec. V. In Sec. VI we summarize the present 
results as well as earlier results for the High cuprates 
in the framework of the f-sum rule. 



II. MODELS 

A. TM model 

Wc first consider a model appropriate for a transition 
metal (compound), referred to as the TM model. Each 
transition metal atom has a five-fold degenerate (n = 5) 
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level. It couples to the other atoms via hopping matrix 
elements t^^, where v = (to, i) is a combined label for a 
orbital index to and a site index i. Thus the electronic 
Hamiltonian is 



^0 / ^ C^fjC^iu 



(2) 



where cj^ creates an electron in the state As discussed 
in the introduction, we consider two different models 
where the atoms are put on a bcc or an A15 lattice, 
describing a transition metal (Nb) or an A15 compound, 
respectively. As discussed above, in the case of the A15 
compound we only consider the transition metal atoms 
and, for instance, neglect Sb in NbaSb.^® This is referred 
to as the Nbg model. 

To describe the hopping integrals, we essentially follow 
Harrison,^*^ and assume that the radial part of the inte- 
grals has a power dependence on the separation of the 
atoms. However, instead of the power five, used by Har- 
rison, we use the power 3.6, more appropriate for Nb.^^ 
Using Harrison notation for the radial part between two 
atomic d energy levels, 



2^1.6 



Vdd,s = Vdd^s 



TO 



R,|3.6. 



,3.6 



(3) 



where 'qdd..c = -16.2, ridd,n = 8.75 and r]dd,s = and to 
is the electron mass. The parameter has been chosen 
in order to reproduce the band with as obtained from 
LDA calculations for Nb^^^, namely rd = 0.7. Since the 
atoms vibrate, their separation can occasionally become 
very small. To avoid that the hopping integrals then be- 
come very large, we have introduced the term containing 
ao in the denominator. We use = 2 A. Eq. (3) shows 
the distance dependence. In addition there are angular 
factors, depending on which TO-quantum numbers are in- 
volved, as described by Harrison. In the model of Nb 
we only consider nearest neighbor hopping, while in the 
A15 model (Nb^) also second nearest neighbor hopping 
is included, since the second nearest neighbors are not 
much further away then the nearest neighbors. 

We consider the case when the phonons couple to the 
hopping integrals (HI). The phonons are approximated 
as Einstein phonons. The frequency Uph = 0.014 eV was 
obtained from the average frequency of Nb metal. For 
each Nb atom we introduce one such phonon in each co- 
ordinate direction. The a;-coordinate of atom i is then 
given by 



2Mu). 



ph 



(4) 



where R^^ is the unperturbed x-coordinate of the atom 
i, creates a phonon in the x-direction on site i and 
M is the mass of a Nb atom. These vibrations couple to 
the hopping matrix elements. 

To obtain the conductivity we calculate the current- 
current correlation function. This requires a definition 



of the matrix elements of the current operator. In our 
model Hamiltonian approach, it is not appropriate to cal- 
culate these as expectation values of the current operator 
between some basis functions, since the basis functions 
underlying our model Hamiltonian are not explicitly de- 
fined. Instead one can use charge and current conser- 
vation, i.e., the requirement that the change of density 
inside some small volume is equal to the current entering 
this volume. This leads to the result 



where /U = (to, i) and v = (m 



(5) 



B. Ceo model 

We next consider a model appropriate for alkali-doped 

fullerenes, referred to as the Ceo model. In these systems 
the till band is partly occupied, and we therefore consider 
a model with a three-fold degenerate fi„ orbital on each 
Cgo molecule i. These orbitals are connected by near- 
est neighbor hopping matrix elements. For the electronic 
part we therefore use the same form of the Hamiltonian 
as above (2), but the orbitals are now three- fold degen- 
erate and placed on a fee lattice. 

The hopping integrals arc obtained from a tight- 
binding description. ^^'^'^ For each of the 60 C atoms in 
a C60 molecule we introduce one 2p orbital pointing ra- 
dially out from the molecule. We then generate orbitals 
of tiu character by forming a linear combination of the 
60 2p orbitals. The hopping between the tiu orbitals on 
different molecules is then determined by the hopping be- 
tween 2p orbitals on different molecules. The 2p orbitals 
couple via a and tt hopping integrals. We use 



= -K/4 



(6) 
(7) 



where Vq = 9.85 eV, do = 1.43 A and L = 0.505 A. 

The calculations were performed for the lattice parame- 
ter 14.24 A. In most calculations we take into account^"^'^^ 
the orientational disorder^^ of the Ceo molecules. 

The important electron-phonon coupling is due to the 
intramolecular phonons of Hg symmetry. There are eight 
such phonons in Ceo, each one being a five- fold degener- 
ate Jahn- Teller mode. Here we only include one degen- 
erate Hg mode per site. We use the Hamiltonian 



^el-ph 
+ 



(8) 



7=1 ia m=lm' = l 



where the phonon coordinate for a phonon with 

quantum number 7 on site i, g is an overall cou- 
pling strength and V^^\ are dimensionless coupling 
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constants^^'^^ given by symmetry. The dimcnsionless 
electron-phonon coupling constant is given by 



(9) 



where N[ii) is the density of states per spin, orbital and 
molecule at the Fermi energy. The current matrix ele- 
ments arc given by Eq. (5) with = R°. 

As a comparison, we also consider a Ceo model where 
the intermolecular phonons couple to the hopping inte- 
grals (HI coupling), instead of the LE coupling consid- 
ered above. This coupling is obtained by displacing the 
molecules from their ideal positions of the fee lattice due 
to the excitations of intermolecular phonons. For large 
values of T, the molecules come unrealistically close to 
each others in our semiclassical theory, neglecting the 
strongly repulsive interaction for small separations, and 
the hopping integrals become unrealistically large. For 
this reason we introduce a modification of the hopping 
integrals between the 2p orbitals in the case of the HI 
coupling . The exponent e~('^~'*°)/^ is replaced by 



-(di-cio)/i£ 



{di-dO)/L I p{d2-do)/L 



Q(d-dQi)/L gid2-do)/L 



(10) 



where di = 3.1 A is the separation of the nearest C atoms 
on neighboring molecules in the equilibrium position and 
^2 = 2 A. For d ^ d2, the hopping integrals are essen- 
tially unchanged, and for d — di and they exactly un- 
changed, while for d <C (i2 the hopping integral is cut off 
at a value which is factor 10 larger than in equilibrium. 



III. METHODS 
A. Quantum Monte-Carlo method 

To establish the properties of our models, we use 
a quantum Monte-Carlo (QMC) approach. For these 
models, the QMC method has no so-called sign problem, 
thanks to the absence of a repulsive Coulomb interaction. 
In the calculation of response functions for imaginary 
times there are then only statistical errors which can be 
made arbitrarily small by improving the sampling. These 
response functions are analytically continued to the real 
frequency axis by using a maximum entropy method. "'^^ 
Although it is nontrivial to control the errors in this 
method, it should still be quite accurate for the response 
functions considered here, due to the simple form of their 
spectra. Thus we are able to quite accurately establish 
the large T behavior of the resistivity for models with 
coupling to phonons. 

In the QMC approach used here,^^ the starting point 
is the partition function 



where Tr is a trace over all states. An imaginary time 
T is introduced, < r < /? = l/iksT). The partition 
function can then be expressed as a functional integral 
over the phonon coordinates as a function of r. For given 
values of the phonon coordinates, the electronic part of 
the Hamiltonian is a one-particle Hamiltonian. The elec- 
tronic degrees of freedom can then be integrated out and 
be expressed as a determinant. Finally, the phonon co- 
ordinates are sampled in a Monte Carlo approach. 

For the LE coupling, the phonons are local and only in- 
fluence the levels on the molecule of the phonon. For the 
Ceo model, this corresponds to a 3 x 3 block in the deter- 
minant obtained in the approach above. The change of 
the determinant when one phonon coordinate is changed 
can then easily be obtained in an updating approach. 
For the HI coupling, on the other hand, each phonon 
influences the hopping integrals to the neighbors of the 
atom of the phonon. Different phonons then couple to 
partly "overlapping" blocks. It is then not possible to 
introduce the simple block form used in the Cgo model. 
This leads to a substantially more complicated updating 
approach, which is discussed in appendix A. 



B. Semiclassical method 

While the QMC method above is very useful in estab- 
lishing the properties of our models, its complexity means 
that it is hard to interpret the results. We therefore 
introduce a much simpler method, where the phonons 
are treated semiclassically. We demonstrate that this 
method is quite accurate for the TM model with HI cou- 
pling, by showing that it agrees quite well with the accu- 
rate QMC calculations. For the Ceo model with LE cou- 
pling, the accuracy is less good, in particular for large T. 
The method is, nevertheless, useful for the interpretation. 

We consider a large super cell with L unit cells, K 
atoms per unit cell and a total of A'^ = KL atoms. Pe- 
riodic boundary conditions are used. Each phonon co- 
ordinate is given a random displacement according to a 
Gaussian distribution centered at zero and the width 



where 



nB{T) = 



n 



Mw. 



'ph 



-ns(T) 



Qhuph/{kBT) _ I ' 



(12) 



(13) 



Z = TVe-^/^, 



(11) 



is the occupation of the phonon mode. In this way, a 
set of displaced coordinates are obtained. These define 
a one-particle Hamiltonian for the electrons. In the case 
of HI coupling, we simply calculate the hopping matrix 
elements using the displaced atomic positions. For the 
LE coupling, we insert the phonon displacements in Eq. 
(8) . Since the coupling contains a factor the Hamil- 
tonian is independent of M for a given A and ojph in the 
case of the LE coupling. 
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To calculate optical conductivity, wc find the eigen- 
states \l) and eigenvalues e; of this Hamiltonian. The 
optical conductivity is then given by 



u' 

(14) 



where O is the volume per atom and /; is the Fermi func- 
tion for the energy e; . The prefactor two comes from the 
summation over spin. We have assumed that the system 
is isotropic, so that it is no limitation to consider the 
conductivity in the x-direction. 

0.15^ , 
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FIG. 3. Resistivity p{T) as a function of temperature T for 
Nbs. The figure compares the semiclassical (broken (A'^ — 36) 
and full (iV = 648) curves) and QMC (circles, = 36) calcu- 
lations. The figure also shows the small (Eq. (25)) and large 
(Eq. (21)) temperature results. The figure illustrates that 
the resistivity of the TM model saturates at large T. Com- 
parison with the QMC results, shows that the semiclassical 
calculation is quite accurate, at least for large T. 
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FIG. 4. Resistivity p(T) as a function of temperature T in 
the ordered Ceo model for ujph = 0.00001 eV and A = 0.6. The 
figure compares the QMC (full curve), the semiclassical (dot- 
ted curve) and the Boltzmann (broken curve) results. The 
phonon frequency was chosen to be so small that p{T) ~ T 
in the Boltzmann theory for all T of interest. 

Fig. 4 compares the semiclassical theory (dotted 
curve) with the QMC (full curve) and the Boltzmann 
(broken curve) theories for the Ceo model with LE cou- 
pling, assuming ordered Ceo molecules. The small T 
behavior is discussed in detail in Sec. VF. Here we 
just notice that the semiclassical theory agrees with the 
Boltzmann theory for very small T and that it agrees 
approximately with the QMC results for small and in- 
termediate values of T. There is, however, a qualitative 
disagreement for large T. The reason is that the strong 
static diagonal disorder introduced by the phonons in the 
semiclassical theory for large T leads to localization. This 
is discussed in more detail in Sec. VH. While the semi- 
classical theory for the Ceo model with LE coupling is 
sufficiently accurate to analyze the results for small and 
intermediate values of T, it is less accurate than for the 
TM model with HI coupling, in particular for large T. 
This is further discussed in Sec. VH. 



Fig. 3 compares the QMC (circles) and semiclassical 
(broken curve) methods for Nbg with iV = 36 atoms in 
the super cell. The QMC calculation has been limited to 
rather large values of T, which is the range of particular 
interest here, and which is also the range of T where the 
calculation can be performed with a reasonable numerical 
effort. The figure illustrates that the semiclassical calcu- 
lation is quite accurate at large T for the TM model. By 
comparing the semiclassical calculation for iV = 36 and 
N = 648 we also illustrate that at large T the result does 
not change much if the size of the super cell is increased. 
For small values of T, however, the discreteness of the 
levels for = 36 would prevent a reliable semiclassical 
calculation for this super cell size. 



IV. RESULTS 

A. TM model 

The full curve in Fig. 3 shows the semiclassical re- 
sults for the Nb3 model. It illustrates how the resistiv- 
ity shows a very pronounced saturation already at quite 
small temperatures. The calculated resistivity at large T 
agrees rather well with the experimentally results, e.g., 
about 0.12 mQcm at T = 900 K (0.08 eV)} This agree- 
ment with experiment is important, since, as we discuss 
below, our saturation resistivity (Eq. (21)) essentially 
only depends on the nearest neighbor distance, the or- 
bital degeneracy n, the filling and the lattice structure. 
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This illustrates that our TM model is appropriate for 
describing resistivity saturation. For small T, the resis- 
tivity grows slower than what is found experimentally, 
which is probably due to the electron-phonon interaction 
being somewhat underestimated in our simple model. 

Fig. 5 compares the semiclassical results for Nb with 
experimental results. The figure shows a surprisingly 
good agreement between theory and experiment, given 
the simplicity of the model and the absence of adjustable 
parameters. The figure illustrates that saturation also 
happens for Nb, but at a much larger temperature scale 
than for Nbg. The reason for this difference is discussed 
in Sec. VD. 



0.12 




0.2 0.3 0.4 

T(eV) 

FIG. 5. Resistivity p{T) as a function of temperature T for 
Nb according to a semiclassical calculation. The figure com- 
pares the semiclassical (full curve) calculation for N — 640 
with experimental results (circles). '^'^ It shows the small (Eq. 
(25)) and large (Eq. (21)) temperature results. The figure il- 
lustrates that there is saturation also for Nb at large T in good 
agreement with experiment. The large T result ((Eq. (21)) is 
slightly exceeded for largo T, since this result is approximate 
and becomes weakly T dependent in a careful analysis (Sec. 
VE). 



der of the Cgo molecules. It is interesting that this "T- 
indcpendent" scattering mechanism gives rise to a weak 
T-dcpendence. The reason for this are discussed in Sec. 
VE. 
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FIG. 6. Resistivity p(T) as a function of temperature T 
and electron-phonon coupling A for the Ceo model according 
to QMC calculations. The phonon frequency is LUph = 0.1 eV. 
The figure illustrates the lack of saturation. For A — 0.80 the 
onset of superconductivity can be seen as a sharp downturn 
in p{T) as T is lowered, due to superconducting fluctuations. 
For A — 1.06 and 1.32, the resistivity has a negative slope for 
small T, indicating an insulating system. The x shows the 
resistivity due to orientational disorder. 

The results for RbaCgo in Fig. 2 were measured at a 
constant pressure and show an approximately quadratic 
dependence on T. If these results are converted to a con- 
stant volume measurement, however, an approximately 
linear dependence on T is found down to T ~ 100 — 200 
K. In agreement with this, Fig. 6 p{T) shows a rather 
linear dependence for A < 0.8 until the superconductiv- 
ity fluctuations set in. The reason for this behavior have 
been discussed earlier. 



B. Ceo model 

Fig. 6 shows QMC calculations for the resistivity of 
the Ceo model according to the QMC calculations. It il- 
lustrates that there is no sign of saturation. Actually the 
curves tend to bend slightly upwards for large T. The 
X indicates the resistivity due to the orientational dis- 
order. This T — resistivity was calculated from Erj. 
(14), i.e., independently of the QMC formalism. The 
curve for A = 0.80 shows signs of superconductivity at 
small T, since the curve tiirns sharply downwards as T 
is lowered, due to superconducting fluctuations. For a 
still larger value of A the system becomes an insulator, 
as illustrated by the negative slope of p{T) for small T. 

The solid curve show the result for A = 0. In this case 
the resistivity is entirely due to the orientational disor- 



C. Comparison of HI and LE coupling 

The results for the TM and Ceo models differ drasti- 
cally. While the TM model shows saturation, the Ceo 
model does not. It is interesting to ask to what extent 
this is due to a difference in the electron-phonon coupling 
(HI versus LE coupling) and to what extent it is due to 
other differences, such as the size of the unit cell, the 
lattice structure and the band width. For this reason we 
have also studied the Ceo niodel assuming a HI coupling. 

The HI coupling in Cgo is due to intermolecular 
phonons, describing the rigid vibrations of the Ceo 
molecules relative to each other. The coupling to these 
phonons has usually been assumed to be weak."^" This 
is also what we find here. We therefore artificially in- 
crease the coupling until A becomes the same as for the 
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intramolecular coupling. Since A ^ Lj~f^ for intermolec- 
ular phonons, we can obtain the increased coupling by 
artifically reducing the phonon frequency LOph- Experi- 
mentally, the intermolecular frequencies fall in the range 
from zero and up to almost 7 meV.^^ We have used a 
value of (jjph = 1.8 meV which is substantially smaller 
than the average frequency of the experimental spectrum. 
The resulting A ~ 0.6 should therefore be substantially 
larger than the experimental value. 

We compare the resistivity in semiclassical calculations 
for the Ceo model with LE and HI coupling in Fig. 7. 
The same values of A ^ 0.6 and LUph = 1.8 meV were 
used in both cases. The molecules are orientationally 
ordered. While the resistivity shows now sign of satura- 
tion for the LE coupling (full curve), the model with HI 
coupling shows a weak saturation (broken curve). This 
becomes even more pronounced if we neglect the rather 
trivial temperature dependence of the Fermi-functions in 
Eq. (14). The resistivity then becomes almost constant 
for HI coupling and large T (dotted curve). For the TM 
model we find a change of slope in p(T) for both HI and 
LE coupling, but the change is more pronounced for HI 
coupling. 
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HI Tp-O 






LE , ..^ 
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0.02 0.04 0.06 0.08 0.1 0.12 0.14 



T(eV) 

FIG. 7. Resistivity p(T) as a function of T for the Ceo 
model considering coupling to the level energies (full line, LE 
coupl.) and coupling to the hopping integrals (broken line, 
HI coupling) according to semiclassical calculations. The Ceo 
molecules arc ordered. The figure also shows results for the 
case when the temperature Tp of the Fermi functions in Eq. 
(14) is put otiual to zero. The figure illustrates that there is 
a large difference between LE and HI coupling for the Ceo 
model. 



V. DISCUSSION 

A. Loss of Drude peak 

We mainly focus on temperatures which are so large 
that the Drude peak is essentially lost. The Drude peak is 



related to intraband transitions between states with simi- 
lar k-vectors. In Appendix B we illustrate that for NbsSb 
in the semiclassical approximation, k-conscrvation is lost 
already at rather small values of T and that the con- 
cept of intraband transitions becomes rather ill-defined. 
Indeed, for large values of T, it becomes a good approx- 
imation to assume that all states couple with the same 
strength via the current operator to all other states, as 
is illustrated in Appendix C and in Fig. 8. The Drude 
peak is then completely lost. Fig. 8 shows that for NbsSb 
the Drude peak is almost completely gone at T = 0.1 eV. 




h (o/w ti (o/w 



FIG. 8. The optical conductivity as a function of the fre- 
quency u) for the (a) A15 and (b) Ceo models in the semiclas- 
sical calculation. The frequency has been scaled by the T = 
band width W . (a) also shows (broken curve) the result of 
approximating all current matrix elements by their average 
(Eq. (C4)). 

B. f-sum rule 

In the large T limit, the f-sum rule provides a very use- 
ful tool for analyzing the resistivity. For model Hamil- 
tonians of the type considered here, the f-sum rule takes 
the form^^ (for a derivation, see Appendix D) 

i^^Mi. = -i^(0|TdO), (15) 

where Tr- is the kinetic energy operator, d is the near- 
est neighbor distance and Q is the volume per atom. As 
discussed above, we assume that T is large enough that 
the Drude peak has been smeared out and that a{u) is a 
smooth function. We furthermore assume that cr(tj) = 
for h\Lo\ > W, where W is the band width. This is ex- 
actly true in the semiclassical treatment and approxi- 
mately true in the QMC treatment. If <t{uj) = a{0) for 
?l|a;| < W, the integral on the left hand side of Eq. (15) 
would be Wa{0) and (t(0) would simply be given by this 
integral divided by W. This is shown schematically in 
Fig. 9. For a more general shape of cr(a;) we write 

a(a; = 0) = -J- / a{uj)nduj, (16) 
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where 7 depends on the shape of cf{lu). To estimate 7 
we assume a certain density of states (DOS) N{e) and 
constant matrix elements of the current operator, as dis- 
cussed in Appendix C. In Table I we give the value of 7 
for different shapes of A''(e), namely a constant 

[0, otherwise, 



TABLE II. The quantities 7 (Eq. (16)) and a (Eq. (20)) 
for a scmi-olliptical DOS (Eq. (19)) as a function of the 
fractional filling p. The results arc symmetrical around 
half-filling (p = ().:>). 



p 


0.1 


0.2 


0.3 


0.4 


0.5 


a 


.041 


.070 


.090 


.102 


.106 


1 


2.63 


2.19 


2.02 


1.93 


1.91 


0:7 


.108 


.153 


.182 


.197 


.202 



a Gaussian 



N{e) 



and a semi-elliptical 



-(2e/WY 



(18) 



N{e) 



otherwise. 



(19) 



DOS. The Tabic illustrates that 7 docs not depend 
strongly on the shape of the DOS. In the following we 
assume a semi-elliptical DOS. 



0(0) 



a(0)/Y 




FIG. 9. Schematic picture of cr(a;). The average over the 
band width is given by a(0)/7. 



TABLE I. The quantity 7 (Eq. (16)) and a (Eq. (20)) 
for a constant (Eq. (17)), a Gaussian (Eq. (18)) and a 
semi-elliptical (Eq. (19)) density of states (DOS) and for 
half-filling. 





Constant 


Gaussian 


Semi-elliptical 


a 


0.125 


0.141 


0.106 


7 


1.44 


1.81 


1.91 


a7 


0.180 


0.255 


0.200 



It is also interesting to study the filling dependence. 

This is shown in Tabic II. The dependence is weak 
around half-filling, but 7 becomes larger for a small fill- 
ing. 



C. Large T behavior 

As above, we consider temperatures which are so large 
that the Drudc peak is gone, but wc furthermore assume 
that the temperatures are small compared with the band 
width. This applies, in particular to many transition 
metal compounds, e.g., the A15 compounds. We consider 
nonintcracting electrons, which should be a reasonable 
assumption for broad band transition metal compounds. 
To apply the analysis above, we have to calculate the ki- 
netic energy Tk- Since T <C VF, we can assume T = in 
the calculation of Tr-. We find that 



rK = 2n [ t 

J-W/2 



eN{e)de = -2naWN 



(20) 



is proportional to the band width W and the orbital de- 
generacy. The shape of the DOS N{e) and the filling 
enter via the parameter a. This parameter is given in 
Table I for different shapes of the DOS for half-filling 
and in Table II for diff'erent fillings and a semi-elliptical 
DOS. Inserting the result for Tk in the f-sum rule (Eq. 
(15)) and using Eq. (16), we obtain 



a(0) = 



7ra7 



3 n M 



(21) 



Here waj/S depends on the details of the electronic struc- 
ture and is of the order of 0.2, cJ-^/fi depends on the lattice 
structure (see Table III), but is of the order 1. 

The result (21) is independent of the band width. This 
follows, since the kinetic energy (20) is proportional to 
W and is cancelled by the W in Eq. (16). 

The quantity ne'^ /(hd) has the dimension of a conduc- 
tivity and it contains the essential material parameters 
n and d. For a transition metal compound, with n = 5 
and ^ 3 A, this leads to an upper limit for the resistiv- 
ity of the order of 0.1-0.2 mOcm. This agrees with the 
saturation resistivities observed for these systems. 

TABLE III. The quantity n/d^ for different lattices, where 
d is the nearest neighbor distance and Q is the volume per 
atom. 





fee 


bee 


A15 


sc 




^=0.707 


573=0-770 


1 = 1.333 


1 
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Using the definition of the mean free path I in the in- 
troduction (Eqs. (1)), we can convert the conductivity 
in Eq. (21) to a mean free path 



/ = cn^d, 



(22) 



where for simpHcity we have assumed that there is only 
one spherical Fermi surface. For a scmi-clhptical DOS 
and half-filhng c = 0.74 (fee), 0.72 (bee) and 0.60 (A15). 
Thus the quantity en's is close to unity for n = 5, as 
appropriate here. This provides a quantum-mechanical 
derivation of the loffe-Regel condition for weakly corre- 
lated systems. 

In particular for the A15 lattice, the second nearest 
neighbor hopping plays a rather important role. The 
separation (0.612a) is not much larger than for the near- 
est neighbors (0.5a), but there are eight second nearest 
neighbors but just two nearest neighbors. For this reason, 
wc also define a distance d which is a weighted average 
of these distances. As weight factors we use the hopping 
matrix elements. Thus we define 



id'} 



(23) 



where d^^ is the distance between the atoms with the 
orbitals p and n. At T = this increases d from 0.5a to 
about 0.57a for the A15 structure. For a semi-elliptical 
DOS and filling 0.4, this leads to a larger saturation con- 
ductivity and a smaller resistivity of about 0.11 mficm 
instead of 0.14 mQcm if the nearest neighbor separation 
is used. This is in better agreement with the calculated 
resistivity. 

It is interesting to study the filling dependence, indi- 
cated by Table II. We consider Sc, which is the first 
element in the 3d series. According to a band structure 
calculation, Sc has about 1.8 3d electrons. Compared 
with a system close to half-filling, such as NbsSb, we then 
expect the saturation resistivity to be about a factor of 
1.5 larger. If we take into account the second nearest 
neighbor hopping, the geometrical factor < d"^ > /f2 is 
similar for the A15 compounds and Sc, suggesting that 
the filling dependence is the dominating factor. Indeed, 
while the saturation resistivity is estimated to be 0.15 
mflcm for NbaSb,^ it is well over 0.2 mficm for Sc,'^^ in 
agreement with the expectations. Similar results arc also 
found for Y.^^ For the other end members of the 3d, Ad 
and 5d series clear saturation does not seem to have been 
observed. 

In a similar way we can use the f-sum rule to estimate 
the resistivity for the Cro model, although the assump- 
tion T <C is now much more questionable, as discussed 
in Sec. VE. Considering a fee lattice, using 7 = 1.91 and 

= 10 A, we obtain 



piT) 



0.288 



Tk{T)/{NW) 



mficm. 



(24) 



Using the band width W = 0.6 eV and obtaining Tk{T) 
from semiclassical calculations for the Ceo model, we 



find the saturation resistivity 0.4 mficm. The calcu- 
lated A = and T = resistivity (0.29 milcm) is below 
this value, while the results for larger values of A and T 
strongly exceed the saturation resistivity. The reasons 
for this are discussed in Sec. VE. 



D. Small T behavior 



In view of the discussion above, we expect the resistiv- 
ity to have an upper limit for models with noninteracting 
electrons scattered by phonons, unless T is very large. In 
many metals, however, the resistivity increases so slowly 
with T, that the corresponding conductivity is much large 
than the limit (21) even at the melting temperature. The 
issue of whether or not the resistivity saturates is then 
not raised. It is therefore of interest to study the low T 
behavior of p{T). For T larger than some fraction of ujph 



we expect 



18 



p{T) = Sn- 



XTks 



(25) 



where fcs is the Boltzmann constant and A is the dimen- 

sionless clcctron-phonon coupling constant. For the TM 
model with HI coupling we define A = A(/x, /x), where 



A(e, e 



dH 

dRia 



(26) 

\^5{ei - e)S{s'i - e - Uph), 



where K is the number of atoms in a unit cell, the a 
summation is over the three coordinates, M is the atomic 
mass, \l) is an eigenstate of H and i labels the atoms in 
the unit cell, flpi is the plasma frequency 



{hnpif 



—y 

37r2 ^ 



^/k[^M^n^-EF). (27) 



where e„k is the energy of a state with the band index n 
and the wave vector k and Ep is the Fermi energy. Clpi 
depends on the average Fermi velocity. 

The straight lines corresponding to Eq. (25) and Eq. 
(21) are shown in Fig. 3. If these lines cross in the 
experimentally accessible temperature range we expect 
saturation. 

It is now interesting to compare our models for Nb and 
Nb^. We obtain similar values of A for the two cases, 
A = 1.0 (Nb^) and A = 0.9 (Nb). A larger value of 
A = 1.7 for Nbg was estimated by AUen'^ while a rather 
similar value was obtained for Nb (A = 1.0) from ab initio 
calculations.^^ We observe that A ~ l/w^^^ depends quite 
sensitively on LOph. Since wc have replaced the whole 
phonon spectrum by three Einstein phonons per atom, 
obtained as the average of the phonon spectrum of Nb,^^ 
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one should not expect very accurate values of A in our cal- 
culation. For the plasma frequency we obtain Qpi = 3.6 
eV (Nbg) and 8.2 eV (Nb), in rather good agreement with 
ab initio calculations 3.4 eV (Nb^)^^ and 9.5 eV (Nb).^^ 

The difference in values of flpi for Nb| and Nb alone 
then leads to a difference by a factor of five in the slope 
of the line from Eq. (25). As a result NbJ shows a very 
pronounced saturation already at small T, while Nb only 
shows sign of saturation at rather large T. The differ- 
ence is due to the fact that Nb3 has a large unit cell with 
many bands and many forbidden crossings. This leads 
to quite flat bands and to small electron velocities. The 
result is a small plasma frequency (Eq. (27)) and a steep 
line from Eq. (25). 

An even more dramatic example is a-Mn, which has a 
unit cell with 58 atoms. One should therefore expect a 
very small plasma frequency and a correspondingly early 
saturation. Indeed, it is found that the resistivity satu- 
rates at about T = 60 K.^^ 

In view of the discussion above, Fig. 8a and Eq. (25), 
it is tempting to write 



a{uj = 0, T) 



Sn'^XTkB 



(28) 



where the first term describes the Drude peak (Eq. (25)) 
and the second term is the conductivity in Eq. (21) at 
saturation. This formula is correct for small T and for 
T which are so large that the Drude peak is gone but 
very much smaller than the band width. Eq. (28) is the 
"parallel resistor" formula of Wiesmann et al..^° 



E. Very large T behavior 

We have so far discussed temperatures which are so 
large that the Drude peak have been washed out, but 
which are small compared with the band width. We now 
focus on values of T which are large enough that the cou- 
pling to the phonons causes a substantial change in the 
band width. Such effects are not very important for typ- 
ical transition metal compounds, which have large band 
widths. They are, however, of substantial interest for the 
Ceo model, for which the fluctuations in the level posi- 
tion become comparable to the band width at values of 
T which can be reached experimentally. 

At such large values of T, there is a rather trivial T 
dependence due to the electron temperature, Tp, enter- 
ing in the Fermi- functions of Eq. (14). This can be 
seen by considering the resistivity due to static disorder. 
Although this scattering mechanism is T-independent, 
the resistivity is, nevertheless, T-dependent. Expanding 
the Fermi functions in Eq. (14) in 1/T, we obtain that 
cr(0) ~ 1/T and p{T) ~ T for very large T. A similar 
dependence also enters for the the electron-phonon scat- 
tering, and it tends to mask some interesting differences 
between level energy (LE) and hopping integral (HI) cou- 
plings. In the following, we therefore freeze the electron 



temperature, Tp = 0, and consider the limit of a very 
large phonon temperature, Tg, i.e., we consider a large 
T but replace the Fermi function by © functions in Eq. 
(14). 

The band width entering Eq. (16) can be approxi- 
mately expressed in terms of the second moment of the 
density of states. The same is also approximately true 
for the kinetic energy. We therefore focus on the second 
moment, 

/oo 
N{e)e^ds, (29) 
-OO 

which can expressed in terms of the Hamiltonian 

where N is the number of atoms in the system. 

We first consider the case of the HI coupling. In our 
semiclassical formalism we can write 



SRia{T) 



dRioc 

+ J E E «f _ 5IU^{T)5Rim + (31) 



2 ^ dRiadRjB 

where the summation over i extends over all atoms. Since 
the displacements SRia are random, we can assume that 

{6R^a{T)) - {6R,o.{T)5Rj{3{T)) = {R'), (32) 

where (i?^) = ksT / {Mujlf^). We then obtain 

E HlAT) - E HlAT = 0) = {R') 5: r ^"^^ • 



dRi. 



(33) 



Explicit calculations for the TM model show that the sec- 
ond term tend to partially cancel the first term, while for 
the Ceo model it adds to the first term. As a crude ap- 
proximation we neglect the second term. The first term 
can be approximately related to the electron-phonon cou- 
pling A. Integrating Eq. (26) we obtain 



NW'^X: 



E(^)' 



(34) 



where we have assumed that A(£, £ ) = A. Assuming 
N{ij) = 1/W, we obtain 

Y.^Ut) -Y.^Ut = 0) = nNXW{T = Q)kBT. 

(35) 
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Assuming that (e^) = W'^/12, as is appropriate for a 
constant density of states, we obtain 



W{T) = W{0)^ 



W{T = 0) ' 



(36) 



where chi = 12. Comparison with expHcit calculations 
for the Ceo model shows that a more realistic value is 
CHI ~ 15. 

We next consider the kinetic energy, Tk- As discussed 
in Sec. VB (Eq. (20)), the kinetic energy is closely 
related to the band width via the quantity a. As T is in- 
creased, however, the shape of N{e) changes somewhat, 
and there is not a perfect proportionality between W{T) 
and Tfc{T). This is illustrated in Fig. 10, where the 
curves describing the T dependence of these two quan- 
tities differ slightly. Nevertheless, from Eqs. (15, 16), 
it follows that the T dependence of these two quantities 
largely cancel in the calculation of a{uj ~ 0) and p{T). 
This is illustrated in Fig. 10, where p{T) has only a 
weak T dependence, once the resistivity has "saturated" 
(at about T = 0.06 cV). The remaining T dependence is 
due to the T dependence of a and 7. 



1.5 



p(T)/p(0.06) 
W(T)/W(0.06) 
Tk(T)/Tk(0.06) 
Y(0.06)/Y(T) 




0.5 





0.02 0.04 0.06 0.08 0.1 0.12 0.14 

T(eV) 

FIG. 10. p{T), Tk{T), W{T) and 7(T) divided by their 
T = 0.06 eV values for the case of a a couphng to the 
hopping integrals (HI) in the Ceo model for A = 0.6 and 
i^ph = 0.0018 eV. The figure illustrates that that the T de- 
pendence of Tk{T) and W{T) arc similar and therefore to a 
substantial extent cancel in Eqs. (15, 16) for HI coupling in 
the Ceo model, leading to a weak T dependence of p{T) for 
large T. The electron temperature Tp =0 and only the boson 
temperature Tb is varied. 

We next consider the case of the LE coupling for the 
Ceo model. In this case the second moment is the sum of 
one contribution from the hopping (off-site elements) of 
H^i, in Eq. (30) and one contribution from the fluctua- 
tions of the level energies (on-site terms) in Eq. (30). To 
obtain the fluctuations in the on-site terms, we rewrite 
the interaction terms as 



^el— phon ^ ^ 



imm Timer 



(37) 



where v^^^i is a boson operator. The average of the 
fluctuation in the level position can then be written as 



X 



N{p) 



T 



(38) 



in the limit of a large T. Combining this with the off-site 
term gives 



3A 

N{0) 



T, 



(39) 



where n = 3. Assuming a constant N^e), we estimate 

that TxiT = 0)/N « -2.6{s'^)].%. For a large T, the 
coupling to the phonons leads to large separations of the 
levels, and we can use perturbation theory for calculating 
the kinetic energy. 



T^(r) = 2^^-^ 



1 Nn{e^)T=o 

2 (e^ - e^) 



(40) 



where we have replaced the denominator by an average 
denominator (e^ — e^) and the limitations on the sums 
to occupied and unoccupied states introduce a factor of 
1/4. A simple estimate of {e^ — e^) is obtained by as- 
suming that the levels have the energies ±Ae/2. Then 
the separation of the levels is Ae = 2{e'^)}/^ , where only 
the on-site contribution to (e^) should be included. At 
large T, however, the on-site contribution dominates and 
we have dropped this restriction. Then 



Tk{T) « -2.6iV 



(41) 



where we have used the same prcfactor 2.6 as below Eq. 
(39). This gives a better agreement with the numeri- 
cal results than the prefactor (3/4) derived from the ar- 
guments above, which is substantially too small, as one 
would expect. The averaging in Eq. (40) greatly favors 
small values of the denominator, while our simple esti- 
mate focuses on large values. The estimate in Eq. (41) 
is also a good estimate for T = 0, as shown above, and 
actually for the whole temperature range. As usual, we 
relate the band width to the second moment. Assuming 
a constant DOS, Eqs. (15, 16) give 



a(0) = 



2.67r7 e^d} {e^)T=o 
6v/l2 {e^W ' 



(42) 



where one factor {s'^)^^ comes from the band width and 
one factor from the kinetic energy. Since (£^)t grows 
with T (Eq. (39)), both the kinetic energy and the band 
width work together to reduce a{0) and to increase p{T) 
as T is increased. Thus we obtain 
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p{T) = -^(1 + cleX^) mOcm. (43) 

From the derivation we obtain cle = 12. A better fit 
to the data is obtained from cle = 16. In addition we 
observe that there is also an appreciable T dependence in 
7(T). These results are illustrated in Fig. 11. In particu- 
lar, we notice that W{T) and Tk{T) have the opposite T 
dependence, and therefore work together in the expres- 
sions in in Eqs. (15, 16). This is in strong contrast to 
the case of HI coupling, where the two T dependencies 
largely cancel each other. 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 



T(eV) 

FIG. 11. p(T), Tk{T), W{T) and 7(T) divided by their 
T = 0.06 eV values for the case of a a coupling to the level 
energies (LE) in the Ceo model with A = 0.6 and tOph = 0.0018 
eV. The figure illustrates that that the T dependence of 
Tk{T) and W{T) are the opposite and therefore work in the 
same direction in Eqs. (15, 16) for the LE coupling in the 
C(io model, leading to a strong T dependence of p(T). The 
electron temperature Tp = Q and only the boson temperature 
Tb is varied. 



F. Lack of saturation in the Ceo model 

By using the f-sum rule, we showed in Sec. VB 
that one should expect the resistivity of the alkali-doped 
hdlcrenes to saturate at about 0.4 mficm. Actually, this 
value is almost reach already at T = (0.3 mfi) due to 
the orientational disorder. One can therefore consider 
the Ceo model as a case where saturation has already 
happened at T = 0. 

This can be further illustrated by considering the re- 
sistivity for a model where all the Ceo molecules have the 
same orientation, i.e., a system without disorder. The re- 
sults arc compared with the resistivity expected from the 
Boltzmann equation in Fig. 12. The phonon frequency 
has been chosen very small, so that the Boltzmann equa- 
tion gives a linear behavior for all T of interest. For small 
values of T the Boltzmann equation and the semiclassi- 
cal theory agree. However, when p{T) becomes of the 



order of 0.3 mficm, shortly before saturation might have 
been expected, the two curves start to deviate. At this 
point we may consider the system has having saturated, 
and the theory in Sec. VE of very large T applies. This 
theory also predicts a linear behavior, but not necessar- 
ily with the same slope as at small T. Simple arguments 
suggest that the two slopes might be of the same order 
of magnitude, as found in Fig. 12. The small T slope 
is, however, related to the properties around the Fermi 
energy, while the very large T slope refers to properties 
integrated over all states. The two slopes should there- 
fore not be expected to be the same. The Boltzmann 
equation is not qualitatively wrong for large T in this 
case, but the relatively good agreement for large T is 
somewhat accidental. 



For the disordered Ceo model, the disorder itself leads 
to a resistivity comparable to the '"saturation resistivity, 
and the "very large T" limit in Sec. VE applies already 
for any finite T . This theory predicts that p{T) has a 
linear dependence on T , as is also approximately seen 
(see Fig. 6). The resistivity could be considered to have 
"saturated" , but this concept is meaningless for the Ceo 
model, since the resistivity grows linearly, with a large 
slope, also after "saturation" . 



We observe that the boson character of the phonons is 
important for the arguments in this section and in Sec. 
V E.^^ Because of this, the number of phonons grow with- 
out limit as T is increased, leading to the corresponding 
growth in the phonon amplitude {R^)- This leads to a 
continuing growth of the band width and reduction of the 
kinetic energy for the case of LE coupling. As a result 
the resistivity does not saturate. 



This is different from the case of electron-electron scat- 
tering, where Fermi occupation numbers enter the the- 
ory. As a result, we have found that there is saturation 
of the resistivity in a simple one-band, symmetric, half- 
filled Hubbard model, at least in the dynamical mean- 
field theory. "'^■^ In view of this, it is interesting that the 
High Tc cuprates are usually considered as examples of 
systems where the resistivity does not saturate, although 
electron-electron scattering is often believed to be the 
dominating mechanism. This issue is addressed in the 
next section. 
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FIG. 12. Resistivity of the ordered Ceo model for 
ojph = 0.00001 eV, A = 0.6 and LE coupling. The figure illus- 
trates that the semiclassical (full curve) and the Boltzmann 
(dashed curve) agree well for small T, but deviates when p{T) 
becomes comparable to the hypothetical saturation resistivity 
(horizontal dotted curve). The figure shows that there is no 
real saturation in this case. 



G. Saturation for High Tc cuprates 

The resistivity in some of the High T,, cuprates is sub- 
stantially larger than one would expect from the loffe- 
Regel criterion. It has therefore been assumed that 
these compounds arc examples of systems where the re- 
sistivity does not saturate. Using the f-sum rule, how- 
ever, we have found that the resistivity saturation is to 
be expected at much higher resistivities than predicted 
by the loffe-Regel criterion or what is found for, e.g., the 
A15 compounds. "'^^ The reason is that the kinetic energy 
is strongly reduced in these systems. This is partly due to 
the strong Coulomb interaction reducing hopping, in par- 
ticular for systems with a small doping x. Furthermore, 
only the — orbital is believed to play an essential 
role, leading to a small degeneracy n = 1. As a result, 
for La2-xSrj:Cu04 we find^* 



P{T) 



0.4 



x{l — x) 



mficm. 



(44) 



This result is much larger than the saturation resistiv- 
ity of the order of 0.1 mficm for the A15 compounds, 
in particular for small x. Experimental resistivities are 
smaller than Eq. (44) , but for small values of x not much 
smaller.^ For these cases signs of saturation are indeed 
seen.^ We therefore conclude that the data are consis- 
tent with saturation. Actually, the data show signs of 
saturation when the experimental resistivity comes close 
to the expected saturation resistivity (44). 



H. Relation to Mott's minimum conductivity 

Within the semiclassical theory, the phonons cause a 
static disorder. The problem discussed here therefore has 
some relations to the conduction in disordered system. 
Thus the LE and HI couplings correspond to diagonal 
and off-diagonal disorder, respectively. While the dis- 
ordered systems are usually studied for small T, we are 
here interested in the large T behavior. In the semiclassi- 
cal theory, however, apart from causing disorder, T only 
enters via the Fermi-functions, and it does not play an 
important role for the qualitative behavior. Below we 
therefore compare our work with the treatment of disor- 
dered systems. 

Diagonal disorder can lead to an Anderson metal- 
insulator transition at T = 0."*^ For the case of off- 
diagonal disorder, however, Antoniou and Economou^^ 
have found that there is no metal insulator transition if 
the Fermi energy is located in some finite region around 
the middle of the band. Our semiclassical calculations 
agree with these results, i.e., we find localization for LE 
but not for HI coupling as T is increased. 

In the QMC calculation of the resistivity, however, we 
see no sign of localization for LE coupling, just a lack of 
saturation. This is natural. Localization depends sensi- 
tively on the phase factors, which are not destroyed in the 
elastic scattering in an disordered system. In the inelas- 
tic scattering by phonons at finite T these phase factors 
are, however, lost, and localization is not expected.*^ The 
effects of the inelastic scattering is properly included in 
the QMC but neglected in the semiclassical treatment, 
and therfore localization shows up in the semiclassical 
but not in the QMC treatment. 

Mott''^ has argued that as the disorder increases, there 
is a discontinuous transition from a metal to an insulator 
at T = 0. He therefore introduced the concept of the 
minimum conductivity 



e2 

CTmin = 0.03—, 

nd 



(45) 



where d is the nearest neighbor atomic distance. Later 
work has argued that the transition from a metal to an 
insulator actually is continuous, but that fTmin may still 
have some relevance for low but nonzero temperatures.*^ 
We therefore make a comparison of (Tmin to the resistiv- 
ity in the TM and Ceo models. Converting Eq. (45) to 
a resistivity, we obtain 



1.6d mficm, 



(46) 



where d is measured in A. Based on experiment, Mott 
deduced a somewhat larger minimum conductivity for 
systems containing transition metal atoms, resulting in 
the maximum resistivity 



1 mOcm. 



(47) 
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Mott derived his result for diagonal disorder. His re- 
sult can most naturally be compared with our saturation 
resistivity for HI coupling (off-diagonal disorder), since 
saturation is most pronounced in this case. The resis- 
tivity pmax is much larger than the saturation resistivity 
obtained above (Eq. (21)) for the TM model with a five- 
fold degenerate orbital (n = 5). For a fee lattice and a 
half-filled semi-clliptical band it takes the form 

Psat = mficm, (48) 

n 

which is of the order of 0.1 mOcm. The corresponding 
conductivity is substantially larger than Mott's minimum 
conductivity. 

I. Alternative explanations 

Cote and Mcisel^ proposed an interesting explanation 
of saturation. They argued that the electrons would not 
see phonons with a wave length A that is much longer 
than the mean free path. They therefore assumed that 
an electron can only be scattered by a phonon if Z > A. 
As T is increased and I is reduced, an increasing fraction 
of the phonons become inefficient as scattering sources. 
The result is that p{T) increases much slower than T at 
large T, in rather good agreement with experiment.^ We 
are now in the position to test this assumption. 

Above, we have studied a model with three local Ein- 
stein phonons on each atom, describing the vibrations 
in the three coordinate directions. This is equivalent to 
study Einstein phonons in q-space. We then write the 
displacement of the atom at the unperturbed position 
as 

where j = 1,...,N labels the A'^ q- vectors and a labels 
the three modes for each q- vector. The corresponding 
phonon amplitude is Ujq. We perform a calculation 
where the phonons are treated semiclassically as before, 
but where the amplitudes treated as random vari- 

ables. This gives the same resistivity as before. We then 
gradually turn off the long wave length phonons, putting 
the corresponding amplitudes Uja = 0. For small T we 
expect this to reduce the resistivity. For large T, how- 
ever, the arguments of Cote and Meisel^ suggests that 
this should not influence the resistivity if A > Z for the 
phonons turned off. 

We group the q-vectors with equal length in shells. 
Shells with q-vcctors of similar length are further 
grouped together in such way that each group contains a 
similar number of q-vectors. Then the groups of phonons 
are successively turned off. The results are shown in Fig. 
13. The figure illustrates that as a group of phonons is 
turned off there is a drop in the resistivity. This is not 



only true for small T but for all T studied here. Con- 
sider for instance the curve with all phonons included 
and T = 0.4 cV. The resistivity p ~ 0.1 milcm corre- 
sponds to Z ~ 3.5 A. The theory of Cote and Meisel then 
assumes that all phonons with A > 3.5 A can be turned 
off without p changing. The figure illustrates that this 
is far from the result of our calculation. This illustrates 
that also phonons with a relatively long wave length con- 
tribute substantially to the large T resistivity, although 
A> L 

Fig. 13 illustrates that phonons with a very long wave 

length make a small contribution to the resistivity for any 
T. The reason is that a long wave length phonon does not 
change the relative separation of two neighboring atoms 
very much, which means that the corresponding hopping 
matrix element is not changed very much. 
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FIG. 13. The resistivity of Nb as a function of T. The 
scattering from phonons of successively shorter and shorter 
wave lengths A is suppressed. For the uppermost curve all 
phonons are considered. In the lower curves the phonons cor- 
responding to the M shortest q-vectors (longest wave lengths) 
were suppressed, where M is marked at the curve. The fig- 
ure illustrates that the long wave length phonons contribute 
about equally much both to the small T resistivity and large 
T resistivity. 

It has also been argued'*^ that resistivity saturation 
can be understood in a Holstein model, somewhat simi- 
lar to our Ceo model. For small T and large A the Hol- 
stein model shows an "excess" resistivity. Similar effects 
are observed in our Ceo model, as is seen in Fig. 6 for 
A = 0.80. The result is that the slope of the p{T) curve is 
reduced as T is increased. To analyze this, we compare 
the calculated p{T) with the resistivity 

p{T) = 0.29 + 17Ar mncm, (50) 

in Fig. 14. The value 0.29 comes from the orientational 
disorder and the term ~ AT is the type of behavior we ex- 
pect for a normal nonsaturating system (e.g., from Boltz- 
mann theory). The slope was adjusted to the results for 
A = 0.26. For such a small value of A there is no sign of 
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saturation in Fig. 6. If the system shows saturation for 
larger values of A, we would then expect the calculated 
resistivity to be below Eq. (50). We find, however, that 
QMC results for large T stay above these results for all 
values of A that we have studied. In the figure this is 
illustrated for A = 0.8. As pointed out in Ref. 46, the 
resistivity in this model actually does not saturate, and 
it was concluded that "saturation" is a misnomer. As 
we have shown above, however, the TM model is a much 
better model of saturation, both because it is much more 
realistic for systems showing saturation, and because it 
also gives results much more similar to experiment. 

In a semiclassical treatment of the type used by Millis 
et al. the "excess" resistivity for large A and small T 
is due to the formation of a highly anharmonic potential 
well for the phonons. This leads to a larger vibration am- 
plitude and an increased resistivity. Similar results are 
found in our QMC calculation, as discussed above. In a 
more realistic model, the electrons would couple to many 
phonon modes, each typically with a substantially weaker 
coupling. Even if the total A may be large, each phonon 
would in such a model have a more harmonic potential 
well, and we would not expect a large "excess" resistivity. 
This further supports our belief that this type of model 
is not appropriate for describing resistivity saturation. 
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FIG. 14. Resistivity p{T) as a function of temperature T 
and electron-phonon coupling A for the Ceo model according 
to QMC calculations. The phonon frequency is ujph = 0.1 eV. 
The straight lines show the resistivity p(T) = 0.29 + ITAT 
mfJcm, wiiere 0.29 is the resistivity due to the orientational 
disorder. The figure illustrates that there is some "excess" 
resistivity at moderate T and large A but no saturation for 
this model. 



VI. SUMMARY 

We have studied models of weakly correlated transi- 
tion metal compound (TM model) and of alkali-doped 
fuUerenes (Ceo model) . These models were studied using 
Quantum Monte-Carlo (QMC) and semiclassical meth- 



ods. The results, as well as earlier results for the High 
Tc cuprates, were analyzed by using the f-sum rule. We 
assumed that T is so large that the Drude peak has been 
smeared out. Then (Eqs. (15,16)) an approximate lower 
limit to cr(0) is given by 



P{T) 



= a{0) 



1 



{u))du) 



\Tk{T)\ 
dW{T) 



(51) 



where Tk{T) is the kinetic energy, W{T) is the band 
width and d is the nearest neighbor distance. 

We first considered T <^W. For the TM model of 
noninteracting electrons, it then followed that Tk ^ W. 
This leads to the simple upper limit 



fid 



(52) 



for the resistivity, where n = 5 is the orbital degeneracy 
of the d-level. This agrees rather well with the satura- 
tion resistivity of many transition metal compounds, and 
it corresponds to a mean free path I ^ d. 

For the High Tc compounds, the kinetic energy is 
strongly reduced by correlation effects. There is a strong 
reduction in the hopping probability of a hole to a neigh- 
boring site if there already is a hole on this site. This 
leads to \Tk\ ^ x(l — x), where x is the doping. The 
corresponding upper limit for the resistivity is then 



he 



e^x(l — x) 



(53) 



where c is the distance between two Cu02 planes. Since 
essentially only the — orbital is involved, the degen- 
eracy factor is n = 1. This resistivity is much larger than 
for the TM model, both because of n = 1 and because of 
factor x{l — x). This limit is therefore apparently never 
exceeded for any high- Tc compound. There are only a few 
cases where the resistivity gets close to this limit, and in 
these cases the resistivity shows signs of saturation. 

Whether or not saturation is actually observed, de- 
pends on how rapidly the resistivity grows for small T's. 
In this limit we have p{T) AT/f^^^ for the TM model. 
For the A15 compounds, e.g., NbaSn, A is fairly large and 
flpi is very small, due to the large unit cell and the quite 
flat bands. The result is that the resistivity grows very 
rapidly for small T and gets close to the limiting value for 
rather small T. The resistivity then shows a pronounced 
saturation. For Nb, on the other hand, ^Ipi is much larger 
and the resistivity grows much more slowly with T, and 
there is only a weak saturation. For most metals, the 
limiting resistivity would only be reached far above the 
melting temperature, due to the slow increase of p{T) for 
small T. 

We also considered very large values of T, where T be- 
comes comparable to the band width. Then both Tk and 
W have strong T dependences. It is important to distin- 
guish between the case when the phonons couple to the 
level positions (LE coupling) and to the hopping integrals 
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(HI coupling). In the former case, Tk decreases with 
T, since the different levels have different energies, and 
hopping is reduced. In the latter case, Tk is increased, 
since the square of the hopping integrals increases with 
T. In both cases W increases with T. In the LE case, 
both effects work together (Eq. (51)) to reduce (t(0) and 
to increase the resistivity. In the HI case, on the other 
hand, the two effects partly compensate each other, and 
the increase in the resistivity is smaller. 

These considerations are very relevant for the Ceo case. 
Due to the orientational disorder, the saturation limit can 
be considered to have been reached already for T = 0. 
Because of the the small band width, however, the T 
dependence of the band width and the kinetic energy be- 
come very important. Furthermore, the coupling is of the 
LE type, so that the T dependence of these two quanti- 
ties cooperate in increasing the resistivity. The result is a 
drastic increase in the resistivity, beyond the "saturation 
resistivity" , and little or no sign of saturation. We may 
therefore consider Ceo to belong to a different class than 
the A15 and High Tc compounds. 
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FIG. 15. Resistivity of Cu, NbsSb (multiplied by a fac- 
tor 1/5),^ Lai.93Sro.o7Cu04 (multiplied by 1/100),® and al- 
kali-doped Ceo (multiplied by 1/100) .'^ The figure also shows 
our estimated saturation resistivities for the latter three cases, 
with the Ceo saturation resistivity (~ 0.4 mf2cm) barely vis- 
ible at about ~400/100=4 /nficm. The figure illustrates that 
the resistivity saturates for NbsSb and Lai.g3Sro.o7Cu04 but 
not for alkali-doped Ceo- 

This is illustrated in Fig. 15, which shows the resis- 
tivity for Cu, NbaSb, Lai.g3Sro.o7Cu04 and hole-doped 
Ceo, where the resistivities of the latter three metals 
have been reduced by factors 5, 100 and 100. The re- 
sistivities of NbsSb and Lai.93Sro.o7Cu04 stay below the 
expected saturation resistivities, while the resistivity of 
Ceo is far above the "saturation" resistivity, shown in the 
lower left corner of the figure. This suggests that the sys- 
tems studied here fall in three different classes, namely i) 
weakly correlated transition metal compounds, showing 
saturation in agreement with the loffe-Regel condition, ii) 



strongly correlated high- Tc cuprates, showing saturation 
but at much larger values than predicted by the loffe- 
Rcgel condition, and iii) alkali-doped Ceo compounds, 
showing no saturation. 
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APPENDIX A: QMC FOR THE TM MODEL 



The Hamiltonian of the TM model can be written as 



ri=l 



(Al) 



where ry labels a given ordering of the Ni, = zNn^ bonds, 
z being the number of atoms connected to a given site 
by the Hamiltonian operator. 

Defining At = 13/ using Trotter decomposition at 
lowest order and breaking up the Hamiltonian in 
terms, the partition function is,^'' 



(A2) 



Integrating out the electron degrees of freedom^^ leads 
to 



Z=[dei{l+BLBL-i...Bi)f 



with 



Bi= ri ri ^ 



-AtH„ 



(A3) 



(A4) 



The matrices 6jj have dimension Nn and have the follow- 



ing form: 
/I 







\0 



.. cosh(Arff:^) 
.. sinh(-Arff^) 






sinh(-ATif^) 
cosh(ATif,,) 







(A5) 



It can be shown^^ that the electron Green function is 
written as: 



g={l + BL...B2Bi)-' 



(A6) 



During the simulation, g and g ^ are constantly stored 
and updated. 
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A Quantum Monte Carlo move is a displacement of a 
phonon coordinate for a given slice. The move is then ac- 
cepted or rejected according to the Metropolis algorithm 
which involves the calculation of the square determinant 
ratio between the electron Green functions after and be- 
fore the displacement, i?^ = [det(5'')/ det(g)] . 

Without loss of generality let us suppose that an atom 
i is displaced in the first slice (so that we can omit 
the higher index in 6^). This will involve a change in 
Bi B[ or Nc = zn^ changes in the factors: 



The first step is to define the new Green function gj 



bl=K^i. k = l,2,...,Ne 



(A7) 



with {ii < ?2 < < ijVe}- In the case only one bi^ factor 
is changed {Nc = 1) the determinant ratio can be easily 
obtained as: 



R = 



det{l + BL...B2bN,:.bi,^h-bi) 



det{l + BL...B2Bi) 
= det[l + (l-5i)(A,,-l)] 



(A8) 
(A9) 



where gi = (1 + bi-^...biBL...B2bpf^...bi-^) ^ is a modified 
electron Green Function and is obtained from g as: 



91 = {bii-i...bi)g{bi^-i...bi) ^ 



(AlO) 



The matrix (Aj — 1) is symmetric and has only four 
matrix elements different from zero, as can be seen from 
Eq. (A5,A7), so that the products in Eq. (A9) can be 
performed in order Nn operation. 

So far it is known how to calculate the determinant as 
long as a single bond is changed. In the more complicate 
case of several bonds, the problem can be reduced to this 
simpler one by noting that the determinant is expressed 
as: 



R = RNc,Nc-iRNc-1,Nc-2---Ri,0 



(All) 



and Rj,j-i is the ratio between two determinants having 
changed only the first j and j — 1 bonds respectively, 



R 



det(l + Bt...B2b,^:,...b,^A,^...b,^A,^...b,■) 
det{l + BL...B2bN,...bij_,Ai._,...bi,Ai,...bi) 

(A12) 



Each of these Nc determinant ratios is given by Eq. (A9) 
with the Green function g replaced by the new one 

= (l + 6i,_i...6^,_^...6^^...6;B^...B26i...6i,)-i 

(A13) 

which has only the first j — 1 bonds updated. 

Once the determinant Rjj-i has been obtained, it is 
necessary to update the Green function Qj-i to the new 
one 'gj which will be used to evaluate Rj+i,j ■ This update 
is done in two steps and requires the knowledge of gj^ so 

that the function g^^ has to be bookkeeped during the 
simulation. 



as: 



gj = (l + 6i,_i...6^^_^...6^^...6ii?i...i?26L...6;^.)-^ (A14) 

gj differs from gj^i only by the substitution bi. b[,. 
It can be obtained using the Green function updating in 
the simpler case of a single bond change^^, namely 



~g,^[9j\ + i9j\-l){X^-i)]-' 



(A15) 



The matrix A = {gjli — l){A.i. — 1) is zero everywhere 
a part from two columns. As a consequence, Eq. (A15) 
can be efficiently performed with the Shermann-Morrison 
formula^^ applied to 'gj^ so that the calculation of gj in- 
volves order {Nnf' operations. 

The second step is then to obtain from gj the Green 
function gj as follows: 

9j = (6i,+,-i...6i,+i6^,)ft(6i,+,-i...6i,+i6^.)-^ (A16) 

Once 'gj\s known it is clearly possible to obtain Rj+ij 
following the same steps we have outlined before. 

For a given Trotter slice and a given phonon coordinate 
the algorithm can be summarized as follows: 

1. Displace coordinate i?j R[ and identify the 
bonds ii < i2 < ■■■ < iN^ which will be affected by 
the atomic displacement. 

2. Set go = 9 and g^^ = g~^, compute gQ and ^ 
using Eq. (A16) and the similar one for g^^. 

3. Perform loop j = 1, ...,Nc over the previously iden- 
tified bonds. 

4. Calculate the matrix Aj^ 

5. Calculate Rjj-i using gj-i and Eq. (A9). 



g, ^ using Eqs. 



6. Update gj_i gj and gj^-^ 
(A15,A16). 

7. End loop over j. 



8. Compute R and check if the proposal move is ac- 
cepted. 



9. If the proposal is accepted update 5jVc-i 
from Eq. (A15). 



9n, 



After the proposed displacement for atom i has been 
accepted by the Metropolis condition, the most straight- 
forward way to proceed would be to obtain the new Green 
function g' (eq. A6), with all the bjj factors updated, as 



g' = {bi„...bi) gN^biM^-bi) 



(A17) 



and then from step 2 of the algorithm obtain the new 
g'g for the atom j = i -\- 1. Note anyway that these two 
steps can be efficiently condensed in one if a particular 
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order for the sites is chosen. If the sites are ordered in 
such a way that ii increase monotonicaUy with i, e.g. 
{li < 2i < ... < TVci}, then Eq. (A17) becomes: 

~9o = {h^^_-K)~^~9NAh^,-K) (A18) 

involving 2ji^ products by 6,, factors less than the most 
straightforward procedure. 



APPENDIX B: LOSS OF MOMENTUM 
CONSERVATION 

At large T the phonon vibrations become very large. 
In the semiclassical treatment of the phonons, this tends 
to destroy the periodicity and therefore it tends to vi- 
olate momcntiim conservation within the electronic sys- 
tem. Below we test how this violation increases with T in 
the TM model using a HI coupling. Qualitatively similar 
results are, however, obtained also in the other models. 
We first calculate the states of the Hamiltonian at T = 0. 
The system is then perfectly periodic and all the states 
\nk,T = 0) can be labelled by a wave vector k and a 
band index n. We use a unit cell with six Nb atoms and 
the band index therefore runs over 30 states. Next the 
states at a finite T are calculated. These states \l, T) are 
labelled by an index I. These states can be expanded in 
the complete set of T = states 

\l, T) = J2 \nK T = 0)(nk, T = 0\l, T), (Bl) 

kn 

For a given state we determine the amount of k-character 
cu{l)=J2\{nKT = 0\l,T)f. (B2) 

n 

or the amount of mixing with states having the band 
index n 

k 

Prom normalization it follows that J^k^k^ — 1 
En'^n =1. We define 

Afe(0=nkE[4'¥, (B4) 

k 

and 

A„(0=30 5][cW]^ (B5) 

n 

where is the number of allowed k-vcctors and 30 is 
the number of band index. If a weight of a given state 
I is equally distributed over nu:/m different k- vectors, 
Afe(Z) = m. In particular, if all effects of periodicity 
are lost, we expect that Ak{l) = 1, since we then ex- 
pect all rik k-components to have equal weight (m = 1). 



On the other hand, if a state contains only one k-vector, 
Afe(/) = rik- Typically in the periodic system, several 
states with different k- vectors are degenerate, e.g., states 
with k and -k may be degenerate. Even at a very small 
amount of disorder, a state of the disordered system is 
then typically a linear combination of states with several 
different k- vectors, and Ak{l) is reduced correspondingly. 
We consider a super cell with periodic boundary condi- 
tions. The value of rik then depends on the size of the 
super cell. For a given amount of disorder, we expect 
that a given state will contain k-vectors from a certain 
fraction of the Brillouin zone. The number of k-vcctors 
increases with the size of the super cell. However, m in- 
troduced above should stay roughly constant. Thus we 
find that definition (B4) gives results which are rather in- 
dependent of the super cell size for values of T which are 
not very small. On the other hand, for T « 0, this defi- 
nition gives results which grow roughly linearly with rik- 
The definition is, however, sensible for the range of T of 
interest here. In a similar way it follows that A„(?) = 1 
if the conservation of the band indices is completely lost. 
We average over all states 

Ai = ^ Aj (1) /{Nn) i = kovn. (B6) 

I 

Fig. 16 shows Aj, for Nbj^ and Nb, where A is an average 
over A(i). The line A = 1, corresponding to a complete 
loss of periodicity, is also shown. The figure illustrates 
that for Nb^ much of the periodicity is lost already for 
T ~ 200 - 300 K. For Nb this happens at higher T, but 
also in this case periodicity is lost fairly quickly. 

The rapid loss of periodicity for Nbg can be related to 
the many flat bands. This means that there are states 
with all k-values within a rather small energy range. 
Then only a small perturbation is needed to mix all these 
different k-values, implying a loss of momentum conser- 
vation. 

In a similar way. Fig. 17 shows that the meaning of 
the band indices is lost relatively quickly for Nug as T 
is increased. This means that the meaning of intraband 
and interband transitions start to loose their meaning. 
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FIG. 16. The quantity A^^ in Eq. (B6) for Nbg and Nb as 
a function of T for nk = 256 allowed k-vectors. Ak measures 

the loss of periodicity. The horizontal line {Ak = 1) repre- 
sents complete loss of periodicity. The figure illustrates the 
rapid loss of periodicity and momentum conservation for the 
Nba model, while this loss happens more slowly for Nb. 



is a Lorentzian. The function j{e, £ ) is shown in Fig. 18 
for two values of T, using the broadening 7 = 0.01 eV. 
The figure illustrates that the function j{e,e ) has only 
a moderate dependence on the energies for T = 0.043 
eV=500 K. We therefore now work out the consequences 
of assuming that the matrix elements of the current can 
be replaced by their average. 
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FIG. 17. The quantity A„ in Eq. (B6) for Nbg as a func- 
tion of T for rik = 256 allowed k-vectors. A„ measures how 
the conservation of the band index is lost, with the horizontal 
line (A„ = 1) showing a complete loss. The figure illustrates 
how the meaning of the band indices is lost relatively rapidly 
for the Nbs model. 



APPENDIX C: CONSTANT CURRENT MATRIX 
ELEMENTS 

In view of the rapid loss of momentum conservation, 
illustrated in Appendix B, it is interesting to consider 
the limit where momentum conservation is completely 
lossed due to the disorder. This is the opposite limit to 
the traditional Bloch-Boltzmann treatment, where the 
scattering is assumed to be so small that k is a useful 
quantum number. In the complete disorder limit studied 
hero, all states arc coupled to all states via the current 
operator. The calculations for the Nbg model show that 
these assumptions, taken literally, are not satisfied. We 
note, however, that the expression in Eq. (14) for the 
optical conductivity can be rewritten as 



j deN{e) J de 



x[f{e)-f{e')]5{ru^-e +e), 



N{e )j{e,s) 



(CI) 



where 



(C2) 

N{s) is the density of states per atom, orbital and spin 
and n is the orbital degeneracy. L{s) = {-y / it) / {e^ + 7^) 
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FIG. 18. Average j{e,e ) of the current matrix elements 
over states with similar energies (Eq. (02)) for Nh^. The 
units are arbitrary. The figure illustrates that for T = 0.0043 
eV=50 K j{e,e ) varies strongly with the energies while for 
T = 0.043 eV=500 K this variation is much less pronounced. 



This average is defined as 
1 



Jav 



{Nnf 



(C3) 



where \i) are the Nn eigenstates of the Hamiltonian. The 
expression (14) for the optical conductivity can then be 
written as 



occ unocc 



—6{hu) ■ 



ei,+ei), (C4) 



where we have considered w > and assumed that 
T << so that we can replace the Fermi functions by 
©-functions. Fig. 8a compares the actually calculated 
a{uj) with the result of (C4), assuming a semi-cUiptical 
DOS (Eq. (C4)). The good agreement for large T gives 
further justification for the assumptions behind Eq. (C4) . 
This gives 



0) 



2TrNrfh 



Jav 



(C5) 



where fj, is the chemical potential. We then need to find 
a relation between jay and N{n), which is obtained from 
charge and current conservation. We first rewrite jav as 



Jav 



1 



(Nn)- 



■ElHi.l/^)!' 



(C6) 
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where |^) is a basis state in a local representation. We 
then use the charge and current conservation in Eq. (5), 
relating the current and hopping matrix elements. This 
gives 



E 

a=x,y,z 



IHialM)!' 



-t 



and for an isotropic system 

.2 _ 1 e2rf2 

Jav 



2 



(C7) 



(C8) 



To relate jl^ to iV(e), we introduce the second moment 

/oo 
N{e)e^de, (C9) 
-oo 

where N(e) is normalized to unity. This quantity can be 
related to the hopping integrals 



(CIO) 



We assume a specific form for N{e), calculate {e"^) for this 
form and then relate it to N{n). Table IV shows results 
for different shapes of the DOS. The table illustrates that 
there is not a drastic dependence on the shape of N{e). 
In the following, we focus the scnii-clliptical DOS, which 
is probably the most realistic one of the three cases con- 
sidered. 

TABLE IV. The quantity {e^}N{^if for a constant (Eq. 
(17)), a Gaussian (Eq. (18)) and a semi-elliptical (Eq. (19)) 
density of states (DOS) and for half-filling. 





Constant 


Gaussian 


Semi-elliptical 




= 0.083 


^=0.159 


^ = 0.101 



Expressing Y^i'^^ix in terms of (e^) in Eq. (C8), we can 



rewrite Eq. (C5) as 



(Cll) 



where Vt/d^ is shown in Table III. The quantity /{hd) 
has the unit of conductivity and Eq. (Cll) can be rewrit- 
ten as 



1 



a(0) 



n/d^ d ^ 



(C12) 



where d is now expressed in A. As seen in Tables IV and 
III, {e'^)N{nf ~ 0.1 and Q./d^ ~ 1. For a transition 
metal, we may use d ^ 3 A and n = 5, which leads to 
p ~ 100 |J^^lcva.. Such a resistivity is indeed typical for the 
saturation resistivity of a transition metal compound. 



APPENDIX D: DERIVATION OF THE f-SUM 
RULE 

In this appendix we derive the f-sum rule, essentially 
following Maldague.^2 We introduce the position opera- 
tor 



The current operator can then be expressed as 



(Dl) 



3x = —[iix,tl\- (D2) 
For w > 0, the optical conductivity is written as 



En — Eq 



where \n) is a many-body state with the energy £"„. By 
inserting Eq. (D2) in one of the two matrix elements of 
jx, one obtains 

-/ a{uj)d^ = -—^m[H,R.],R.m- m 

Performing the commutators, we find 

Y}[H, Ra],Ra] = E d%t^^■^dai^^.a, (D5) 



where d,^^ is the distance between the sites with the or- 
bitals v and /x and a labels the coordinate. This result 
is true for noninteracting systems as well as interacting 
systems of certain types, e.g., with an on-site Hubbard 
interaction. We now assume only nearest neighbor hop- 
ping, replacing di^^ by d. Furthermore, we assume the 
system to be isotropic, so that all directions a are equiv- 
alent. For a three-dimensional system, the commutator 
on the right hand side of Eq. (D4) is then one third of 
the result in Eq. (D5). This gives 



a{u)du) • 



(D6) 



where Tk is the kinetic energy. For a two-dimensional 
system the factor 3 in the denominator is replaced by a 
factor 2. This result can also be generalized to a finite 
temperature. In the case of the TM model, however, the 
atomic separations cannot be treated as constants, since 
they vary as the phonons are excited. The coordinates 
in Eq. (D4) can then not be taken outside the average 
(...). We can, nevertheless, recover an expression like Eq. 
(D6) by defining an appropriate average separation d{T). 
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